Slow stress relaxation in randomly disordered nematic elastomers and gels 
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Randomly disordered (poly domain) liquid crystalline elastomers align under stress. We study 
the dynamics of stress relaxation before, during and after the Polydomain-Monodomain (P-M) 
transition. The results for different materials show the universal ultra-slow logarithmic behaviour, 
especially pronounced in the region of the transition. The data is approximated very well by an 
equation a{t) ^ aeqis) -\- A/{1 -\-a\og t). We propose a theoretical model based on the concept of 
cooperative mechanical resistance for the re- orient at ion of each domain, attempting to follow the 
soft-deformation pathway. The exact model solution can be approximated by compact analytical 
expressions valid at short and at long times of relaxation, with two model parameters determined 
from the data. 

PACS numbers: 64.70.Md, 75.10.Nr, 83.80.D 



Liquid crystalline ordering in a confined, randomly 
quenched geometry has been the subject of consider- 
able research in recent years. Non-aligned polymer- 
stabilised or dispersed liquid crystals and, in partic- 
ular, the nematics in aerogels or other porous media 
have been studied with an eye on the universal 
properties of systems with weak random fields. The ba- 
sic scaling observation, that at long length scales the 
static random fluctuations dominate over the dynamic 
thermal fluctuations, has been verified by rigorous the- 
oretical analysis and experimental studies of correla- 
tions and susceptibilities in a wide variety of systems 

. The critical slowing down of all relaxation processes 
in systems with a random, glass-like order is also a well- 
established universal phenomenon H,^, with a char- 
acteristic stretched-exponential or power-law behaviour 
as opposed to a normal exponential relaxation of over- 
damped ordered systems. Accordingly, it has not been 
surprising to find slow relaxation modes in randomly 
confined nematic liquid crystals However, al- 

though the concept of a random orientational effect im- 
posed by the large amount of inner surfaces seems in- 
tuitively correct, the direct application of continuum 
random-field models to systems with strong anchoring 
on sub- micron length scales is more difficult. 

Nematic elastomers and gels have been a focus of ex- 
tensive recent research for a number of different, but 
equally compelling reasons. In these systems the liq- 
uid crystal ordering and director are coupled to the 
mechanical degrees of freedom - stresses and strains 
of the underlying polymer network. As a result, local 
elastic torques are unbalanced and the non-symmetric 
Cosserat-like elasticity leads to a number of unique 
physical effects such as soft elasticity and mechanically- 
driven orientational switching (see the review for de- 
tails). It has been argued that the network crosslinks 
could act as local sources of quenched disorder and the 
nematic gel is analogous to the spin glass with ran- 
dom magnetic anisotropy jio). The application of ex- 
ternal stress to (polydomain) nematic networks disor- 
dered on large scales results in a critical transition into 
the aligned state with an increasing degree of long-range 
order. This polydomain-monodomain (P-M) transition, 
analogous to the spin-glass alignment by magnetic field. 



has been studied experimentally in some detail [ pJj , 
confirming that the alignment proceeds via the reori- 
entation of correlated regions ('domains') rather than 
the growth of the favoured ones. Because the polymer 
chains within each domain are anisotropic, the rotat- 
ing director causes shape changes, Fig.l, and thus the 
materials can accommodate the external deformation 
without significant stress response. In this Letter we 
study the dynamics of this transition by increasing the 
extensional strain e in controlled small steps and mon- 
itoring the relaxation of stress cr(t). 




FIG.l Schematic illustration of a stretched polydomain elas- 
tomer. As each correlated domain rotates its director towards 
the extension axis, its shape is changing. The shear strain associ- 
ated with the rotation (indicated by shading of a single domain) 
needs to be compatible with the surrounding elastic medium. 

We find this relaxation to be extremely slow, increas- 
ingly so in the vicinity of critical stress for the P-M 
transition, and following the logarithmic law over sev- 
eral decades of time. Similar logarithmic decays has 
been observed in other systems, as different as nemat- 
ics in silica gels and avalanches at the angle of re- 
pose and are reminiscent of the 1// noise kinetics 
|]l3| . We propose a theoretical model for a rate con- 
stant vanishing with an essential singularity due to the 
cooperative mechanical barriers for each domain's ro- 
tation. The resulting kinetic equation, having a form 
0=—m0^e~^^^^ gives a solution which at short times 
resembles the power-law and at long times can be 
approximately interpolated as 1/(1 + 2 log t/t*), with 
the crossover time = (rnv?)~^. Both regimes corre- 
spond well to the experimental results. We expect this 
behaviour to be universal across the randomly disor- 
dered systems where, as in nematic elastomers and gels. 
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the relevant order parameter is coupled to the elastic 
modes: the mechanical compatibility requires the coop- 
erative nature of elastic barriers. 

The materials used in this experiment are three 
different types of side-chain liquid crystalline poly- 
mers: polysiloxanes [Q crosslinked via rather flexible 
tri- functional groups, polyacrylates crosslinked by bi- 
functional chemical groups Q and by 7-radiation [|6| . 
All materials had slightly different crosslinking density, 
around 5-10%. The backbone anisotropy is very differ- 
ent for siloxane and aery late chains, the anisotropy of 
chain radii of gyration being R\\/R±_ 1.6 and 1.06, 
respectively (for Gaussian polymers {R^\ ^) = li\\^±L 
with L the chain contour length). The glass transition 
temperature for polysiloxane elastomer was around -5C 
(with nematic-isotropic transition T^i ~ 42C), while the 
polyacrylates become glasses at Tg ~ 50C (with T^^ ~ 
HOC) [Q. Hence the experiment was carried at 30C 
for polysiloxane elastomer and at around 90C for both 
polyacrylates. In spite of all these differences, the re- 
sults obtained are distinctly similar and in this Let- 
ter we concentrate mostly on the polysiloxane system, 
which has more dramatic quantitative effects due to the 
higher chain anisotropy. The experimental approach is 
straightforward. Strips of polydomain nematic rubber 
(~ 15x5x0. 3 mm) were suspended on a stress gauge and 
extended in controlled strain fashion in a box which was 
thermostatically controlled. We applied a consecutive 
fixed- step extensions of ~0.5 mm every 24 hours (thus 
providing an effective strain rate of s ~3-10~^s~^). The 
data for variation of response force with time has been 
collected and then converted to the nominal stress a 
(and the extension - to engineering strain s = AL/Lq). 
The question of temperature control required a seri- 
ous attention. The accurate measurement of very small 
changes in stress over large time intervals demanded 
the experimental error to be brough down to 0.05% in 
stress and 0.1° in temperature. The partially and fully 
aligned nematic elastomers are very thermally sensitive. 
We, therefore, included an additional correction to ac- 
count for this effect. 

The stress-strain dependence of polysiloxane nematic 
rubber, going through the P-M transition is shown in 
Fig. 2. The sets of data points shown are collected 
at fixed intervals of 20s, 100s, 500s and 24 hours af- 
ter the strain increment from the previous-day value. 
One clearly sees the major trend of the P-M transi- 
tion, exhibited at all times: When the threshold stress is 
reached, the correlated nematic domains begin to rotate 
towards the direction of extension. Because the average 
polymer backbone is prolate along the local nematic di- 
rector n, the aligning of average n results in the effective 
lengthening of the whole sample, thus accommodating 
the imposed strain and creating the stress plateau over 
a large interval of deformations. The strain at the end 
of the equilibrium stress plateau is determined by the 
anisotropy of chain shape R\\/R±. However, the ques- 
tion of what is the equilibrium stress at each value of 
s requires a detailed analysis of its kinetics. The in- 
set in Fig. 2 shows the analogous plot for a chemically 
crosslinked polyacrylate, with a much smaller backbone 



anisotropy and the stress plateau region. 
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FIG. 2 Experimental values of the nominal stress in polysiloxane 
nematic rubber, as observed at 20s (0), 100s (A), 500s (0) and 24 
hours (•) after the strain increment. The numbered points corre- 
spond to the data sets used in Fig. 3. The thick dashed line is the 
saturation level of 'apparent power-law'; the solid line shows the 
true equilibrium values (Teqis) the stress would achieve, ijLgiven 
time for inverse-logarithmic decay to saturate [see Eq.(Bj) and 
Fig. 4]. The inset shows the same plot for chemically crosslinked 
polyacrylate. 
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FIG. 3 Six selected relaxation plots cr(t) for points, labelled from 
1 to 6 in Fig. 2, in polysiloxane sample. Note, that the plots for 
points p3 (e = 0.13) and p4 (e = 0.21) are hardly distinguishable, 
even though they are at each end of the stress plateau. The slow 
(logarithmic) relaxation at long times is evident, with the slope 
significantly increasing as the P-M transition is being approached 
from both aides. Thin solid lines show the 'apparent power law' 
model, Eq.(Bj), which fits only the short-time region of each plot. 

Fig. 3 shows the variation of stress with time after 
the strain increment for several selected points in Fig. 2. 
The logarithmic scale contracts the long-time region 
and enhances the effect of decreasing (j{t). One imme- 
diately recognises the fact that, especially for the points 
on the stress plateau {i.e. during the P-M transition), 
the equilibrium stress has not been achieved even after 
the 24-hour interval. Fig. 4 gives an example of data 
analysis. The relaxation of stress at relatively short 
times accurately follows a power-law with a seemingly 
universal exponent This variation has been repro- 
duced in all three different materials (studied at dif- 
ferent temperatures), for all points on the stress-strain 
curve of Fig. 2. At longer times a logarithmic decay was 
evident. Again, this has been a feature for all materials 
and at all strains, however, the slope of the logarithmic 
decay was very different, dramatically increasing for the 
materials with higher chain anisotropy and towards the 
P-M transition stress plateau. In fact, the logarithmic 
decay was hard to detect for the first and last points in 
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Fig. 2 for highly anisotropic polysiloxane nematic rub- 
ber and, equally, for the chemically crosslinked poly- 
acrylate (which is known to have a very small backbone 
anisotropy). Finally, the elastomer samples at large de- 
formations (of 10-30% in our study) always develop a 
distinct necking near the rigid clamps (which do not 
allow the reduction of sample width required by incom- 
pressibility). As a result, even at deformations long 
past the stress plateau and the samples well-ordered in 
the middle, there are still polydomain regions near the 
clamps are in an incomplete state of alignment (this is 
clearly seen by the eye and under the microscope). Sim- 
ilarly, well before the P-M transition is reached, there 
are internal processes associated with domain wall lo- 
calisation [0 , also contributing to slow cooperative re- 
laxation. Therefore, it is attractive to attribute the 
residual logarithmic decay away from the P-M transi- 
tion to such artefacts of non-uniformity and conclude 
that the main stress relaxation occurs during the P-M 
transition and is due to the rotation of correlated ne- 
matic domains. 
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FIG. 4 Analysis of relaxation for the point p4 (0) where the stress 
has been normalised for convenience to vary between 1 and 0. The 
'apparent power law' model with the exponent (dots) fits the 
data well until the crossover time t* is reached. A good overall 
fit of the stress a(t) is given by the inverse-logarithm function 
c = CTeq + A/(l + Q: log t) , the solid line. 



then proceeds via the director rotation towards its soft- 
elasticity equilibrium, Acr(t) ~ /i(^||/^± — l)^(^), while 
the dynamics of 6 is controlled by the standard model- 
A equation = —fhO^. Note the cubic force, which 
is the direct consequence of soft elasticity and leads to 
the power-law decay 6 = (2m t + ■^)~^/^, where the 
initial misalignment is the function of strain step. 
This power law is indeed observed at early stages of 
relaxation for all materials and deformations studied. 
However, the long-time logarithmic behaviour requires 
a more delicate analysis. 

So far we have discussed the soft response of an in- 
dividual domain. Surrounded by its neighbours, each 
of which has its own director and a different set of soft 
strains, any given domain will face an elastic barrier for 
its relaxation. In an extreme situation when all neigh- 
bouring domains are already aligned and cannot find 
any soft pathway to accommodate the external strain, 
the shear deformation of the given domain will not be 
allowed by mechanical compatibility and no relaxation 
would occur jl^]. Therefore, for each individual domain 
to relax its local stress, one requires a cooperation of 
its neighbours, providing a 'gap' for the required set of 
strains. We shall approximate this effect by estimating 
the effective rate constant m of escape over the average 
non-soft barrier. Let us first introduce the mean angle 
{0{t))^ the average of misalignment angles of all differ- 
ent domains and thus the measure of non-relaxed part 
of the stress. As {0) decays to zero after each strain 
increment, the total barrier in the system increases as 
(^o) — (^(^)) with the initial condition (^o) proportional 
to the strain step. If we assume that this energy is 
evenly distributed between all non-relaxed domains (the 
number of which is ex (^)), the effective escape rate be- 
comes m ~ exp[— /3i'-^^^|^y^], leading to the dynamical 
equation 



l(^) = _^e-"/W {Of 
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To understand the results, one needs to recall the 
basic concepts of soft elasticity in nematic gels (or 
any other Cosserat-like elastic media with an indepen- 
dent Goldstone degree of freedom); the reader should 
refer to ^ for details. Considering a region of lo- 
cally uniform director (a 'domain'), the rubber-elastic 
energy can be reduced to zero if the director rota- 
tion and the Cauchy strain tensor A are related by 

the equation A = :£(6>)^/^:£(0)~^/^, where the uniaxi- 
ally anisotropic polymer step-length matrix is given by 
% = + (^11 — i±)ninj. For a plane stretching 

and director rotation as in Fig.l, when the deformation 
£ = ^zz — 1 is applied, the elastic response can be re- 
duced if the director rotates by ^ y/2£\\/{i\\ - i^e^l^. 
To comply with the soft-elasticity pathway, this rotation 
requires an associated shear Szx ^ ^(^||/^± — 1) sin 2Q 
[H. Therefore, a given domain subjected to an external 
extension e finds itself out of elastic equilibrium, with a 
director misaligned by d and the effective energy density 
A/ ^ - Ife^ at smah l9, where /i ^ ny^ksT - 

lO^J/m^ is the rubber modulus. The relaxation of stress 



(with a small re-arranging of constants). We expect the 
parameter u to be proportional to the strain step and 
to the overall strain which is a practical measure of 
the P-M transition. Fig. 2. Both m and u should also be 
proportional to the chain anisotropy — 1). Eq. 

(H) integrates for times between zero and t and gives 



t 



mu^ 



1 + ProductLog - (1 - i^)e-i+^/^o 



(2) 



with the direct solution d(i) expressed via the special 
function ProductLog[z], a solution for w oi z ■= we^ 
(we have omitted the (..) notation for the mean angle 
S). The constant in its argument is ~ 1/e at expected 
ii/6>o <C 1. It is straightforward to plot the solution (see 
Fig. 5) and verify that the essential singularity in the 
rate constant at ^ ^ 0, Eq. (|^), results in a critical 
slowing down of the long-time decay. 
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FIG. 5 The log-log plot of scaled solution 6/u against scaled time 
t/t* = mu^t, as suggested by Eq.(ph, circles. The arrow shows the 
positien of crossover time with two approximate regimes, given by 
Eqs.(pl), shown by dots ('apparent power law') and solid line (in- 
verse logarithm). The inset shows the model activation parameter 
u vs. £, determined in three ways described in the text: (i) - 0? 
(ii) - and (iii) - •. 

There are several qualitative conclusions to be drawn 
from the Eqs.(|l|)-(|^). At short times, when is not 
small, the exponential in the decay rate is irrelevant and 
we recover the power-law behaviour ~ This 
'apparent power law' becomes invalid at a crossover 
time t* ^ {mv?)~^ (when 6 ^ u)^ after which the 
inverse-logarithm function gives an excellent interpola- 
tion of the exact solution (0): 
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for 
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1 + 2 log(mii^ t) 



for t>r 



(3) 



The corresponding stress relaxation behaviour, a{t) = 
(Jeq + const • 9{t) closely reproduces the data, presented 
in Figs. 3 and 4. The experimental values for crossover 
time t*, slope and saturation level of 'apparent power 
law' and fit of the long-time logarithmic tail provide 
a sufficient number of independent measurements for 
each data set to accuartely determine the two model 
parameters, m and u. The inset of Fig. 5 shows the 
variation of parameter u for the polysiloxane sample 
as function of position on the stress-strain curve, de- 
termined by three different methods: (i) using the dif- 
ference between cjeq and the 'apparent power law' sat- 
uration being ^ 0.27 (ii) identifying the numerator 
of the long-time inverse-logarithm fit with u and (iii) 
combining the slope of the 'apparent power law' m 
and the crossover time for each plot, u = (mt*)~^/^. 
The plateau value u ^ 0.0023 (and the corresponding 
plateau for the parameter m ^ 700) strongly depend on 
the material and temperature, presumably via the chain 
anisotropy, but are also approximately linear functions 
of the small strain step. 

Considering some experimental ambiguities (such as 
the effect of non-uniform necking) and theoretical sim- 
plifications, we consider the proposed model to provide 
a good qualitative description of the data. All three 
materials show the same type of response, with only 
parameters m and u reflecting the difference. Much re- 
mains to be done to build a full understanding of stress 



relaxation phenomena in nematic elastomers with, it 
appears, the unavoidable effects of quenched disorder. 
One needs to study the effect of temperature, in partic- 
ular, close to the nematic transition point T^i. The 
detailed role of chain anisotropy R\\/Ri_ remains to 
be investigated. Correlating the stress relaxation with 
equilibrium autocorrelation functions, expected to fol- 
low the activated scaling ||] as the spin-glasses and ran- 
dom nematics do, would also be desirable. A refinement 
of the theoretical model is desirable, particularly to ad- 
dress the role of small linear (semi-soft) corrections to 
the cubic force. 
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